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FUNCTIONS THROUGH THE CUNTZ RELATIONS 
AND NEW INTERPOLATION PROBLEMS 
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Abstract. We introduce connections between the Cuntz rela- 
tions and the Hardy space H2 of the open unit disk D. We then 
use them to solve a new kind of multipoint interpolation problem 
in H2, where for instance, only a linear combination of the values 
of a function at given points is preassigned, rather than the values 
at the points themselves. 
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1. Introduction 

One motivation for studying representations of the Cuntz relations 
comes from signal processing, sub-band filters, and their applications 
to wavelets. This falls within a larger context of multiscale problems, 
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see for example |T^. In this work we study the Cuntz relations in 
a different context, and introduce connections between them and the 
Hardy space H2 of the open unit disk D. We prove in particular the 
following results: Let be a finite Blaschke product of degree M, and 
let Ci, . . . , Cm be an orthonormal basis of H2 6H2. A function / be- 
longs to H2 and has norm less or equal to 1 if and only if it can be 
written as 



AI 

/m G H2, are uniquely defined, and are such that 

M 

Ih2 = ll/jllHa- 



(1.1) 

where /i, 
(1.2) 



From now on, we denote by II/II2 the norm of an element of H2. Using 
Leech's factorization theorem (see Section [2] below), we prove that, 
equivalently, / belongs to H2 and has norm less or equal to 1 if and 
only if it can be written as 



(1.3) 
where 

(1.4) 



E 



a 



is a Schur function, i.e. analytic and contractive in D. 



Representation (11. ip allows us to solve various interpolation problems 
in H2 by translating them into tangential interpolation problems at 
one point (in fact at the origin) in H2^. The solution of this latter 
problem, or more generally, of the bitangential interpolation problem 
in Hg^'' is well known. See for instance [H [2]. 

Similarly, the representation (II. 3p allows us to solve various interpola- 
tion problems in H2 by translating them into tangential interpolation 
problems at one point (here too, in fact at the origin), for C^^"'"^-valued 
Schur functions, whose solution is well known. See for instance [12] for 
the general bitangential interpolation problem for matrix-valued Schur 
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functions. 



We now illustrate these points. First note that, for b the Blaschke 
product with zeroes the points ai, . . . , um, (H.ip leads to 



M 



1.5) 



/(a,) = 5^e,(a,)/,(0), i 



M. 



For preassigned values of f{ai), i = 1, . . . , M, this reduces the Nevanlinna- 
Pick interpolation problem for M points in H2 to a tangential inter- 
polation problem at the origin for functions in H2^, whose solution, as 
already mentioned, is well known. The novelty in the present paper is 
by exploiting the above reduction scheme to solve multipoint interpo- 
lation problems in H2. For example, consider the following problem: 



Problem: Given M points oi, . . . , um in D and u = (mi U2 
£ixM -feC, find all / G H2 such that 

M 

(1.6) ^M^/(a^)=7. 
Solution using (II. ip : It follows from (11.51) that 

M AI M 

^uJitti) = ^C^Uiej{ai))fj{0) 
j=i e=i 

M 



um) e 



;i-7) 



with 
(1.8) 



M 



^Uiej{ai), j = l, 



M. 



For preassigned value of the left side of (ll.7p this is a classical tangen- 
tial interpolation problem for C*^-valued functions with entries in the 
Hardy space. Let v = (yi V2 ■ ■ ■ vm) & C^^'^^ ■ Assuming vv* 7^ 
we have that the set of solutions is given by 



f2{z) 
\fM{z)J 



VV 



V V 



vv 



\9m{,z)/ 
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where gi, . . . ,gM G H2 and 

M I |2 

Ell Alii = !^+Eii*ii^' 

e=i e=i 

It follows from (11. ip that a function / G H2 satisfies (II. 6p if and only 
if it can be written as 

/(^) = ^Ee,(z)t;;+(ei(z) e,{z) ■■■ eM{z)) B{z) ^'^'^ 
^ \9m{z)J 

where we have denoted B{z) = (^Im + (z — 1)^). Note that B is an 
elementary Blaschke factor, with zero at the origin. 

Solution using (II. 3^ : In the case of representation 01. 3p . we have 
similarly 

M 

(1.9) f{ae) = 5^e,(a,)ai,(0), £ = 1, . . . ,M. 

i=i 

For preassigned values of f{ae), i = 1, . . . , M, this reduces the Nevanlinna- 
Pick interpolation problem for M points in H2 to a tangential inter- 
polation problem at the origin for matrix- valued Schur functions (II. 4p . 
As in the previous discussion we exploit the above reduction scheme 
to interpolation problem in the Schur class for multipoint interpola- 
tion problems. For example, in the case of the interpolation constraint 
ffL6D . it follows from (OD that 

M M M 

(1.10) 

with Vi, . . . ,vm as in (11.81) . For preassigned value of the left side of 
(ll.lOp this is a classical tangential interpolation problem for C^^~^^- 
valued Schur functions. 



Problems of the form (II. 9p have been studied for M = 2, under the 
name multipoint interpolation problem, in [6]. In that paper, an invo- 
lution (fi of the open unit disk which maps ai into 02 is used. Then one 
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notes that the function 

F{z)- 
satisfies the symmetry 

F{^{z)) = JF{z), where J 



1 

1 



This reduces the interpolation problem in H2 to an interpolation prob- 
lem with symmetries in Hg. Unfortunately this method does not extend 
to the case M > 2. For a related interpolation problem (for Nevanlinna 
functions), see also [H], where the n-th composition of the map f is 
equal to the identity map: (f°"'{z) = z. 

In the present paper, we use a decomposition of elements in H2 as- 
sociated with isometries defined from b and which satisfy the Cuntz 
relation. The representation of Hardy functions, proved in |5], plays 
a major role in the reduction to interpolation problems in the setting 
of Schur functions. To ease the notation, we set the discussion in the 
framework of scalar-valued functions, but the paper itself (as well as 
[3]) is developed for matrix- valued functions. Besides being a key player 
in complex analysis, the Hardy space H2 of the open unit disk plays an 
important role in signal processing and in the theory of linear dynam- 
ical systems. An element / in H2 can be described in (at least) three 
different ways. In terms of: (1) power series, (2) integral conditions, 
or (3), a positive definite kernel. More precisely, in case (1) one sees 
/ as the ^-transform of a discrete signal with finite energy, that is the 
^-transform of a sequence {fn)neNo ^ h- 

00 00 

fiz)=J2Uz-, \\f\\h/^-J2\fn\'<00, 
n=0 n=0 

In case (2), one expresses the norm (in the equivalent way) as 



1 



2lT 



^ = — sup / \f{re^')\'dt< 00, 

^TT rG(0,l) Jo 

and sees / as the transfer function (filter) of a £1-^2 stable linear system. 
See [20] • In case (3) we use the fact that H2 is the reproducing kernel 
Hilbert space with reproducing kernel jz^- From the characterization 
of elements in a reproducing kernel Hilbert space, a function / defined 
in © belongs to H2 if and only if for some M > 0, the kernel 

(1.11) -2 /(')/("-)• 

^ ^ 1- zw* M 
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is positive definite there. The smallest such M is ||/||2- For M = 1, 
rewriting f ll.lip as 

a{z)a{w)* - h{z)h{wy ( \ (^ fi \\ 
Y^^, , with a{z) = (1 -zf{z)) , 

and using Leech's factorization theorem (see next section), it was proved 
in [5] that / admits a (in general not unique) representation of the form 

(Ji{z) 



:i-i2) f{z) 



1 - za2iz) 



where a{z) = J is analytic and contractive in the open unit disk. 

Let now a G ©, and 

baiz) 



1 — za* 

In [To] was proved that the map 

(1-13) TJ{z) = /(b.(^)) 

1 — za* 

is from H2 onto itself and unitary. In the present paper we replace ba 
by an arbitrary finite Blaschke product, and define a counterpart of 
the operator T^. If M = deg b, we now have, instead of the unitary 
map, Ta a set of isometrics Si, ... , Sm in H2, defined as follows: Take 
ei, . . . , Cm be an orthonormal basis of the space H2 6H2. Then, 

(1.14) iSjh){z) = e,iz)h{biz)), heU2, 
with Si, ... , Sm satisfying the Cuntz relations: 

M 

(1.15) ^5,5; = /h„ 

Ih2, if 3 = k, 
0, otherwise. 



:i.i6) s*Sk 



It follows from these relations that every element / G H2 can be written 
in a unique way as (11.11) : 

M 

f{z) = Y,e,iz)fM^)), 
i=i 

where the fj G H2 and satisfy ( 11. 2p 



M 
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We note that Si, ... , Sm in (II -14^ form a finite system of M isome- 
tries with orthogonal ranges in H2, with the sum of the ranges equal 
to all of H2. Thus they define a representation of the Cuntz relations. 
This is a special case of a result of Courtney, Muhly and Schmidt, see 
[HI Theorem 3.3]. We send the reader to [18] for a survey of the rele- 
vant literature and in particular for a discussion of the related papers 
[281 [2Z]- For completeness, we provide a proof, in the matrix- valued 
case, using reproducing kernel spaces techniques (see Section Hj). As 
already mentioned, one motivation for studying representations of the 
Cuntz relations comes from signal processing. Our present application 
of the Cuntz relations to Leech's problem from harmonic analysis is 
entirely new. The immediate relevance to sub-band filters is a careful 
selecting of the Cuntz isometries, one for each frequency sub-band, see 
[26| Chapter 9]. In the case of wavelet applications, the number M 
is the scaling number characterizing the particular family of wavelets 
under discussion. 

We note that relations with the Cuntz relations in the indefinite inner 
product case have been considered in [9], in the setting of de Branges 
Rovnyak spaces; see [TSKTB]. This suggests connections with interpola- 
tion in these spaces (see [21 Section 11], ^3\), which will be considered 
elsewhere. We briefly discuss some of these aspects in Section [3 

Our paper is interdisciplinary, a mix of pure and applied, and we are 
motivated by several prior developments and work by other authors. 
This we discuss in section|21in]and|l]below. For the readers convenience, 
we mention here briefly some of these connections. One motivation 
comes from earlier work [T] on two-sided, and tangential, interpolation 
for matrix functions; see also references [2] through [5], and [Hj. We 
make additional connections also to to interpolation in de Branges- 
Rovnyak spaces [13], to wavelet filters, see e.g., [T7], and to iterated 
function systems, see [18j by Courtney, Muhly, and Schmidt, and [28] 
by Rochberg, to Hardy classes [291 EQ] and to classical harmonic anal- 
ysis; see e.g., [221 ESI El ES] • 

The paper consists of six sections besides the introduction. Sections 2 
and 3 are of a review nature: In the second section we discuss Leech's 
theorem and in the third section we discuss the realization result of [5]. 
In Section 4 we consider, in the matrix-valued set of operators 

which satisfy the Cuntz relations, and were considered earlier in [18] in 
the scalar case. Section 5 is devoted to the proof of the matrix version 
of (11.31) . We use the representation theorem of Section 5 to solve in 
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Section 6 new types of multipoint interpolation problems. Finally we 
outline in the last section how some of the results extend to the case 
of de Branges Rovnyak spaces. 

2. Leech's theorem 

As already mentioned in the introduction, we set the paper in the 
framework of matrix- valued functions. When the Taylor coefficients /„ 
are C'^^'^-valued, one defines a C'^^'^-valued quadratic form by 

[/, /] =• 7^ / {f{e^')rf{^')dt = J2 /n7n. 

The space H*'^'^ consists of the functions for which Tr [/, /] < oo In [5] 
a representation theorem for elements / G 112^"^ such that [/, /] < Iq 
in terms of Schur functions was presented. See Theorem 13.11 below. 
Recall first that a C^^'^-valued function a defined in the open unit disk 
is analytic and contractive in the open unit disk if and only if the kernel 

(2.1) x,(...0 = ^^^« 

1 — zw 

is positive definite in the open unit disk. Such functions are called 
Schur functions, denoted by S^p^i, Given a G yv^'i and a C^^^-valued 
function A analytic in the open unit disk, the kernel 

./Nr^/ A{z)A{wy - B{z)B{wY 

A{z)K,{z,w)A{wy = ^ ' ^ y 

1 — zw 

where B = Aa, is positive in D. Leech's theorem asserts that the 
converse holds: if A and B are respectively C'^^'^-valued and C'^^^- 
valued functions deGned in D and such that the kernel 

(2.2) Aiz)A{wy - B{z)B{wY 

1 — zw* 

is positive definite in D, then there exists a G yp^t such that B = Aa. 
Two proofs of this theorem hold. The first assumes that A and B are 
bounded in the open unit disk, and uses a commutant lifting result of 
M. Rosenblum. See PH] for Rosenblum's result and |30i Example 1, p. 
107], |7| for Leech's theorem. The other proof requires only analyticity 
of A and B in D, and uses tangential interpolation theory for Schur 
functions, together with the normal family theorem. One can extend 
these arguments to functions of bounded type in D, or even further 
weaken these hypothesis. For completeness, we now outline a proof 
of Leech's theorem for continuous functions A and B. We first recall 
the following: Let G N and let Wi,. . . ,Wn G D, ^i,. ■ ■ e 
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and r]i,...,r]N G C'^. The tangential Nevanlinna-Pick interpolation 
problem consists in finding all Schur functions cr G J/'p^i such that 

The fact that the function K„{z, w) defined by fl2.1l) is positive definite 
in D implies that a necessary condition for the tangential Nevanlinna- 
Pick interpolation problem to have a solution is that the N x N Her- 
mitian matrix P (known as the Pick matrix) with ij entry 

(2.3) p =§El^U]lIll 

is non negative. This condition is in fact also sufficient, and there are 
various methods to describe all solutions in terms of a linear fractional 
transformation. See for instance [121 [El 1211 [22]. With this result 
at hand we can outline a proof of Leech's theorem as follows: We 
assume given two functions A and B, respectively C'^^'^-valued and 
C'^^'^-valued, continuous in © and such that the kernel 

A{z)A{w)* - B{z)B{wY 



1 — zw* 

is positive definite there. Consider wi,W2. ■ ■ a countable set of points 
dense in the open unit disk. The Hermitian block matrix with (.j entry 



A{w,)A{wjY - B{wi)B{wj] 

1 — WiW* 



,3 = 1,. ..,N, 



is non negative, and therefore, by the above mentioned result on Nevanlinna- 
Pick interpolation, there exists a Schur function G S^p^i such that 

A{wi)aN{we) = B{we), l=l,...N. 

To conclude the proof, one uses the normal family theorem to find a 
function a G such that A{wi)a{wi) = B{we) for i e N. By 

continuity, this equality extends then to all of D. 

3. A REPRESENTATION OF Hg^^ FUNCTIONS 

Leech's theorem can be used to find a representation of elements of H2 
in terms of Schur functions, as we now recall. See [5lll]. The following 
result is proved in [5]. In the scalar case, it was proved earlier by 
Sarason using different methods. See [351, P- 50], [321 |33l [31] • In the 
discussion we shall find it convenient to partition a G S^(p+i^^i as 

with 0"! being C^^'^-valued and (T2 being C^^-valued. 
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Theorem 3.1. Let H G H2^^. Then the following are equivalent: 

(1) It holds that: 

(3.2) [H,H]<I„ 

(2) The kernel 

(3.3) — ^ - H{z)H{wr 

1 — zw* 

is positive definite in D. 

(3) There is a Schur function a e ^(p+'?)x'? (see so that 

(3.4) H{z)=a^{z){I^~za2{z))-\ 

The key to proof of this theorem is to note that the kernel (13.31) can be 
rewritten in the form fl2.2p with 

A{z) = {lp zH{z)) and B{z) = H{z), 

and apply Leech's theorem: There exists a E ^(p+i^^i as in f l3.ip such 
that A{z)a{z) = B{z), that is, 

aiiz) + zH{z)a2{z) = H{z). 

Equation (13. ip follows. 

We note that an extension of the previous theorem to elements in the 
Arveson space was given in |3, Theorem 10.3, p. 182]. 

4. The Cuntz relations in 
Let 6 be a finite Blaschke product of degree M, and let 

j^{b) = H2 e 6H2. 

It is well known that this space is finite dimensional and i?o-iiivariant, 
where 

Rofiz) = l^^l^^, fe^ib). 

z 

Let 

(/i(z) ... fM{z)) =C{Im-zA)-' 

denote a basis of J^{b), where (C, A) G C^^*^ x C^'^^*^ is an observable 
pair, namely: 

n^^okerCA" = {0}. 

Since the spectrum of A is inside the open unit disk, the series 

00 

P = J2 A^*C*CA^ 

i=0 
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converges and P > 0. The matrix P is the Gram matrix (observabihty 
Gramian in control terminology) of the basis /i, . . . , Jm, and satisfies: 

P=^r(/i(e^*) ... /./(e^*))* (/i(e^*) •.. fAiie'')) dt. 

^TT Jo 

This matrix turns to be identical to the Pick matrix defined in fl2.3p . 
We denote by H2 the Hilbert space of C^-valued functions with entries 
in H2, and with norm 




■n|lH2- 

n=l 



We note that is the reproducing kernel Hilbert space with repro- 
ducing kernel , 

We now introduce the Cuntz relations into this framework. This was 
treated earlier in [18] using different methods. 

Theorem 4.1. Let ei, . . . , cm be an orthonormal basis of Jif{b) and, 
forj = l,...M 

(4.1) {Sjh){z) = ej{z)h{b{z)), heHl 

Then the Sj satisfy the Cuntz relations: 



M 



(4.3) S*S, 



^H^, if j = k, 



1 0, otherwise. 

Proof: We proceed in a number of steps. The proof of the Cuntz 
identities f l4.2p -( l4l3|) is given in STEPS 4 and 5, respectively. 

STEP 1: The set H2(b) of functions of the form 

F{z) = f{h{z)), /GHf, 

with norm 



l-^llH^(b) = llH^. 



is the reproducing kernel Hilbert space with reproducing kernel 



b{z)b(w)* 
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This can be checked directly, but is also a special case of jHl Theorem 
3.1, p. 109]. 

STEP 2: The operator Mg. of multiplication by ej is an isometry from 
H2(6) into Hf. Furthermore, the range of Me^ and M^^. are orthogonal 
for j 7^ k. 



Indeed, let u,v E C^. It holds that 



v*u j = k and m = n, 
0, otherwise. 



We use that multiplication by b is an isometry from H2 into itself. If 
n = m and j = k, the claim is clear. If n = m and j 7^ k, this is just 
the orthogonality of cj and Ck- If n > m, we have 

(ej6"M, ekb"'v)iiP = (ej6"-™M, ekv)iiP = 0, 

since ejb"'~''^u G 6H2 is orthogonal to CkV whose components belong to 
J^{b) = H2 6H2. The case < m is obtained by interchanging the 
role of j and k. 

STEP 3: Let Ci, . . . , cm denote an orthonormal basis of Jif{b). Then, 

(4.4) H^ = ©*i,e,Hf(6). 

Indeed, the reproducing kernel is written in terms of the orthonormal 
basis as (see for instance [HI (6) p. 346], [31] ) 

M 

(4.5) K,{z,w) = J2eA^)ie,Hr- 
Thus 



(4,6) r^= , ,1' '-^W'W 

1 — zw* 1 — b[z)b[w)* 1 — zw* 

with 

_ e,(^)(ej(^;))*Jp 
k,{z,w)- ^_^^^^^^^y ■ 

Equality (14. 6 p expresses the positive definite kernel ^ sum of 

positive definite kernels. The reproducing kernel space associated to kj 
is ej'H.2{b). Therefore, (14. 4p holds as a sum of vector spaces; see [HI p. 
352]. Since, by STEP 2, ejH2(6) is isometrically included into H2, the 



M 
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sum is orthogonal. 

STEP 4: Sj and Sk are isometries, with orthogonal ranges when j ^ k. 

The fact that Sj is an isometry follows from STEPS 1 and 2. Indeed 
the range of Sj is in by STEP 2 and 

WSjhWl,. = \\hmn^,i,) (by STEP 2) 
= \\h\\l. (by STEP 1). 

Furthermore, for f,g& and j ^ k, 

{Sjf,S,g)^. = {M,J{h),M,^g{b))^. = 

by STEP 2. 

STEP 5. It holds that , S.S* = Jh^. 

Indeed, by the properties of multiplication and composition operators 
in reproducing kernel Hilbert spaces, we have that, with pw{z) = 1 — zw* 
and u G C^: 

Pw Pb{w)[z) 



M M 



Thus 

(i:5.5;)(^„.)^X:^-^e,,.)(e,H) 

where we have used (14. 5p and (14. 6p . This ends the proof since the 
closed linear span of the functions — is all of H2 . □ 

pw 

Thus we have the following decomposition result for elements in H2. 
When M = 1, this result reduces to the fact that the operator 
defined in (I1.13P is a unitary map from H2 into itself. 

Theorem 4.2. Let h he a finite Blaschke product of degree M , and 
let ei, . . . , Cm be an orthonormal basis of Jif{b). Then, every element 
H G H2^'' can be written in a unique way as 

M 

(4.7) H{z) = Y,eAz)H,m), 
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where the Hj G Hg^'' and 



M 



(4.8) 



Proof: We define operators 5*1, ... , Sm as in (14. ip . Let H e Hf^*^ and 
^ e C^. It follows from the definition (14.11) of the Sj and from (14.21) 
that 

M 

i=i 

where Hj G Hg^'^ is defined by Hj^ = S*{H^). Taking now into account 
(14.31) we have 

M 
M 



M 

M 



and hence (14. 8 p holds. 



□ 



5. Representation of elements of H?^^ 



We now present a generalization of Theorem 13. II To this end, we gener- 
alize (II. 4p to a partitioning of a matrix- valued functions a G ^i'^'^P+i)^i 
as 



(5.1^ 



a 



with (Til, • • • ; CTiM being C^^'^- valued and o"2 being C^^"^- valued. 
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Theorem 5.1. Let b be a preassigned finite Blaschke product, and let 
ei, . . . , Cm be an orthonormal basis of Jif{b). Let H G 'R^^'^ . Then the 
following are equivalent: 

(1) Condition holds: [H,H] < Ig. 

(2) There exists a G yi^^P+i)^i such that 

(5.2) H{z) = |^^e,(.)ai,(6(.))j (/, - b{z)aMz))r\ 
where a G yiMp+ci)xg is as in flETjl . 



Proof: By Theorem 14. 2[ H subject to (13. 2p can be written in a unique 
way as (14. 7p . and it foUows from (14. 8 p that 



M 



Using Theorem 13.11 with the function 



G 



H2 



\HmJ 



uMpxq 
^ -"-2 ) 



we see that there exists a G yi^-^p+<i)><i (see (15. ip ) such that 

(/g - za2{z))~'^. 



(H,{z)\ /au{z)\ 
H2{z) 



\Hm{z)J 
The resuh follows using (14. 7p 



\(riM{z)J 



□ 



The results in [5] are a special case of a family of interpolation problems 
with relaxed constraints. See [23 [21] ■ We plan in a future publication 
to consider these results in our new extended setting. 

6. New interpolation problems 

We have outlined in the introduction the connections between multi- 
point interpolations and representations ( II. ip and ( II. 3p . We now add 
some details. Interpolation problems whose solutions are outlined in 
the present section will be considered in full details in a future publi- 
cation. 
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The case of (11 .ip : We consider the foUowing problem: Find aU func- 

rP> 
*-2 



tions H G Hn^*^ such that 



M 



(6.1: 



for some pre- assigned matrices ^i,...,^a/ £ C*"^^ and 7 G C^''. To 
solve this problem we use (14. 7p with 



h{z) = 

A basis of J^{b) is given by 

1 z 



z — a 
1 — za* 



M 



M-1 



1 — za* (1 — za 



*\2 ' 



1 — za 



(see for instance pT]). Set 
(6.2) 
Since 



^i^) — [l-za'-^P (l-^a*)2-^P {l~za*)M "-P j ■ 



(6.3) 

and with 



h{a) = h'{a) = --- = h^^'-^\a)=Q, 



\HM{b{z))J 



-"-2 ) 



we have that 



H{a) = ^(a)^(O) 
H'{a) = E'{a)J^{0) 



^(Af-i)(o) = ^(A^-i)(a)^(0). 
Therefore, the interpolation problem (16.11) is equivalent to 

cje{o) = 7, 

with C G C^^^P given by 

M-l 

C=Y.i,E^^\a). 

j=0 
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When CC* > 0, this in turn can be solved using [21 Section 7], or 
directly, as 

J^iz) = C*{CC*)-'l + {hip + l)C*iCC*)-^C) ^{z), 

where 5f G ii2^^^'^. The formula for H follows. We note that 

JT] = 7*(CC*)"S+ [^,^]- 
The case where CC* is not invertible is solved using pseudo- inverses. 

The case of (11.31) : We here assume first that p = q = 1 and 

M 

Z - tti 



za} 

where the are distinct points in D. We have now 

C = (l 1 ■■• l) and ^4 = diag (a^, . . . , a^). 
Note that the pair (C, ^4) is observable. Define 

(6.4) P^^ = T^irv*^ i,j = l,...,M. 

Namely we are in the case (12. 3p with the = 1 and the 'qj = 0. In 
other words, P is the Pick matrix obtained while interpolating all the 
points a£ to the origin. We mention the papers [36] and [23] for a 
related discussion. 

Proposition 6.1. Let u = [ui U2 ■ ■ ■ um) G C^^^ and 7 G C 6e 
pre-assigned. Then the following are equivalent: 

(1) It holds that 

M 

''^^Uih{ai) = 7 and ||/i||h2 < 1- 
e=i 

(2) h is of the form 

C{I - zA)-ip-i/Vi(6(^))(l - h{z)aMz)r\ 

where a = \ '^^\ G is such that 

(6.5) (mP1/2 0) (t(0) = 7, 
where P is defined by (16. 4p . 
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When b{z) 



z-a \M 



1-za* 



) for some M G N and a G D, one obtains a 



different kind of interpolation problem, as we now explain. Rewriting 
O as 

(6.6) Hiz)il, - biz)a2ibiz))) = E{z)a{b{z)), 

where E{z) is given by (16.21) with p = 1, and 



aiiz 



Differentiating (16.61) M — 1 times and taking into account (16.31) . we 
obtain that 



(6.7) 



H'ia) 



( E[a) Opxp\ 



E\a) 



-"pxp 



^2(0) 



\E^^^~'\a) Opxp/ 



This allows to reduce (16.11) to a standard tangential interpolation prob- 
lem for Schur functions. 



7. The case of de Branges Rovnyak spaces 

Let s be a Schur function. The kernel ks{z,w) = is positive 

definite in the open unit disk, and the associated reproducing kerenl 
Hilbert space will be denoted by ^(s). Such spaces were introduced 
and studied in depth by de Branges and Rovnyak in [16J. When s is 
an inner function (and in particular when s is finite Blaschke product), 
we have 

^(s) = H2 e SH2. 

In general, ^{s) is only contractively included in H2. Let moreover b 
be a finite Blaschke product. We have 



1 - sibiz))s{b{w))* _ 1 - s{b{z))s{b{w))* 1 - b{z)b{w) 



(7.1) 



1 — zw* 1 — b(z)b{w)* 1 — zw* 

= f e,(.)e,(,„)-i--^"'<=»^<''<'"»' 



1 - b(z)b(w 



with ei, . . . , Cm an orthogonal basis of J^{b). This decomposition al- 
lows us to define the operators Si, ... , Sm as in (14. ip , so that the 
following holds: 
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Theorem 7.1. The operators Si, ... , Sm are continuous from Jif{s) 
into Jif{s{b)) and satisfy the Cuntz relations: 

M 

(7-2) '^SjS* = Ij^(s{b)), 



(7.3) S*S, 



k 



I^{s), if j = k, 
0, otherwise. 



Proof: We proceed in a number of steps. 

STEP 1: The reproducing kernel Hilbert space ^{s, b) with reproduc- 
ingkernel ^~l^!^lfj-jl^(^yf^^ consists of the functions of the form F{z) = f{h{z)), 
with f G Jif{s) and norm 



1-^11.^(^,6) = \\J WJ^Xs)- 



This follows from a direct computation. 

STEP 2: The formula 

(r,(/c,(,)(-,w;))) (z) = K{z,b{w))ej{wy, weB. 

defines a bounded densely defined operator, which has an extension to 
all of Jif{s{b)), and whose adjoint is Sj. 

This follows from the decomposition (17. ip . 

STEP 3: O holds. 

Indeed 

M \ \ M 

^SjS* j ks(b){-,w) j (2;) = ^ej{z)ksib{z),b{w))ej{wy 
j=i J J j=i 

= k,(j,){z,w), 

and hence, by continuity, equality (I7.2p holds in J^{s{b)). 

STEP 4: For j ^ k we have 

(7.4) Cj^is, b) n Ck^is, b) = {0} 



20 D. ALPAY, P. JORGENSEN, I. LEWKOWICZ, AND I. MARZIANO 



Indeed, let H2(6) be as in STEP 1 in the proof of Theorem 14. 1[ We 
have Jif{s) C H2 and hence 

^{s,b) C H2(6). 

Thus fl7.4p follows from STEP 2 of that same theorem. 

STEP 5: Let M^.. denote the operator of multiplication by Cj. It holds 
that 

(7.5) J^{s{b)) = (BjUM,^^{s,b). 

This follows from the decomposition (IT.ip . which implies that the sum 

M 

i=i 

holds, and from STEP 4, which insures that the sum is direct. 
STEP 5: dH hold. 

Indeed, from STEP 4, the range of the operators M^^ and Mg^ are or- 
thogonal for j ^ k, and Mf,. is an isometry. □ 

Finally, we remark that (17. 2p leads to decompositions of elements of 
the space Jif{s{b)) in terms of elements of the space Jif{s) similar to 
(II. ID : Every element / G J^{s{b)) can be written in a unique way as 

M 

/(^) = 5^e,(z)/,(&(^)), 
i=i 

where /i, . . . , /m € J^{s). Furthermore 

M 

\\f\\jf'{s{b)) = 

Multipoint interpolation problems can be also considered in this set- 
ting, building in particular on the recent work of Ball, Bolotnikov and 
ter Horst [13] on interpolation in de Branges Rovnyak spaces. This will 
be developed in a separate publication. 

Our paper is meant as an interdisciplinary contribution, and it involves 
an approach to filters and to operators having its genesis in many differ- 
ent fields, both within mathematics and within engineering. We hope 
that we have succeeded at least partially in communicating across tra- 
ditional lines of division separating these fields. As a result our listed 
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references included below is likely to be incomplete. We thank in par- 
ticular Professor Paul Muhly for improving our reference list. 

With apologies to Goethe and to Frenchmen: 

Mathematicians are like Frenchmen: whatever you say to them they 
translate into their own language and forthwith it is something entirely 
different 

Johann Wolfgang von Goethe. 
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